Astronomy & Astrophysics manuscript no. aa6475 


© ESO 2008 


February 3, 2008 





o 
o 



> 

o 



Estimation of Galactic model parameters in high latitudes with 

2MASS 

A. Cabrera-Lavers'-^, S. Bilir^, S. Ak^, E. Yaz^, and M. Lopez-Corredoira' 



' Institute de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain 
l ^" ^ ^ GTC Project Office, E-38205 La Laguna, Tenerife, Spain 

, Istanbul University Science Faculty, Department of Astronomy and Space Sciences, 34119, University-Istanbul, Turkey 



^ Received XX; accepted XX 

m , 

, ABSTRACT 
^ ■ 

T-H Context. In general, studies focused on the Milky Way's structure show a range of values when deriving different Galactic parameters such 
as radial scalelengths, vertical scaleheights, or local space densities. Those values are also dependent on the Galactic coordinates under 
consideration for the corresponding analysis, as a direct consequence of be observing a structure (Our Galaxy), that is far from being as smooth 
' and well-behaved as models usually consider. 

I Aims. In this paper, we try to find any dependence of the Galactic structural parameters with the Galactic longitude for either the thin disc or 
I the thick disc of the Milky Way, that would indicate to possible inhomogeneities or asymmetries in those Galactic components. 
Methods. Galactic model parameters have been estimated for a set of 36 high-latitude fields with Two Micron All Sky Survey (2MASS) 
t| . photometry. Possible variations with the Galactic longitude of either the scaleheight and the local space density of these components are 
' explored. 

, Results. Galactic model parameters for the different fields show that, effectively, they are Galactic longitude-dependent. The thick disc 
scaleheight changes from ~800 pc at 150° < / < 230° to ~1050 pc at |/| < 30°. A plausible explanation for this finding might be the effect 
. 5^ of the flare in this Galactic component that changes the scaleheight (h-) with Galactocentric distance (R) following the approximate law: 
h-(R) = (940 ± 20) X [1 - (0.12 ± 0.02)(R - Rq)]. The effect of the flare is more stressed in some lines of sight than in others, producing the 
5_j observed changes in the parameters with the Galactic coordinates used to derive them. 



Key words. Galaxy: general — Galaxy: stellar content — Galaxy: structure — Infrared: stars 



1. Introduction 

Detailed study of star counts in the Milky Way permit us to re- 
cover basic structural parameters, such as its disc scalelength 
and disc scaleheight. Our knowledge of the structure of the 
Milky Way, as inferred from star count data, is about to en- 
ter the next level of precision with the advent of new surveys 
such as 2MASS. This will require the refinement of models of 
Galactic structure to fit the parameters of the basic components 
of Milky Way. 

Researchers have used different methods to determine the 



Galactic model parameters. In Table 1 of iKaraali et al.l (120041) 
we can found an exhaustive summary of the different values 
obtained for the sctructural parameters of the discs and halo of 
the Milky Way. One can see that there is an evolution for the 
numerical values of model parameters. The local space density 



and the scaleheight of the thick disc can be given as an example. 
The evaluations of the thick disc have steadil y moved towards 
shorter scaleheights, fro m 1.45 to 0.65 kpc jGilmore & Reidl 



19831 : IChen et al.L 1200 lb and higher local densities (2-10 per 
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cent). In many studies th e range of val ues f or the parameters is 
large. For example, Che n et al.l (120011) and lSiegel et al.l (l2002l) 
give 6.5-13 and 6-10 per cent, respectively, for the local space 
density for the thick disc. 



Different model parameters revealed for the Galactic disc 
may due to the warp and flare. The disc of the Milky Way is 
far from being radially smooth and uniform. On the contrary, 
it presents strong asymmetries in its overall shape. While the 
warp bends the Galactic plane upwards in the first and second 
Galactic longitude quadrants (0° < Z < 180°) and downwards 
in the third and fourth quadrants (180° </ < 360°), the flare 
changes the scaleheight as a function of radial distance. 
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The warp is present in all Galactic componen ts: d ust 
(Drimm el & Spergelll200l[ lMarshall et al., 200 61). gas (iBurt 
198 8; Drimmel & Spergel', 2001; Nakanishi & Sofue', 12003 
Levin e. Blitz & Heiles, 2006; Voskes & Burton, 2006 ), and 
stars dLopez-Corredoira. Betancort-Riio & Beckmanl 12002 : 



Momany et all 120061) . All these components have the 
same node position and their distributions are asymmet- 
ric. However, the amplitude of the warp seems to de- 
pend slightly on the component one looks at: the dust 
warp seems to be less pronounced than the stellar and 
gaseous warps, that share approximately the same am- 



Northern Sky and one at Cerro Tololo Observatory (CTIO, 
Chile) to complete the Southern counterpart of the survey. 
Observations cover approximately 97 per cent of the Sky with 
simultaneous detections in J (1.25;um), H (1.65 yum), and Ks 
(2.17 yum) bands up to limiting magnitudes of 15.8, 15.1, and 
14.3, respectively. 

In this paper we have used data from the All Sky Release 
of 2MASS, made available to the public on March 2003, 
that includes a po int source catalogue of -470 million stars 



plitude ("Lopez-Corred oira. Betancort-Rijo & Beckmanl 2002t 
[Momany et al., 2006). 

The stellar and gaseous flaring for t he Milky Way 



are also compatible dMomanv et all 120061) . showing that 



(ICutri et al.L l2003h . We extracted those areas at Galactic lati- 
tudes of 45° < ^ < 55° and 60° <b< 70°, where both the dif- 
ferential reliability and completeness of the 2MASS catalogue 
are 0.99, thus the nominal limiting magnitudes of the survey 
can be easily achieved. 



the scalehei ght increases with the Galactocentric radius for 

R > 5 kpc dKent. Dame & Faziol Il99ll; IPrimmel & Soergell 2.1. Stellar density from the red clump population 



ergell 
2002t 



2001"; 'Naravan & Jog^. l2002t iLopez-Corredoira et all 1200 
Momany et al., 2006). The behaviour of this flare in the cen- 
tral discs of spiral galaxies is not so well studied due to in- 
herent difliculties in se parating the several contribution s to the 
observed counts or flux. lLopez-Corredoira et al" (2004 1, for ex- 
ample, found that there is a deficit of stars respect to the predic- 
tions of a pure exponential law in the inner 4 kpc of the MiUcy 
Way, that could be explained as being a flare which displaces 
the stars to higher heights above the plane as we move to the 
Galactic center. 

In this scenario, where on the one hand the mean disc 
(z = 0) can be displaced as much as 2 kpc between the 
location of th e maximum an d the minimum amplitude of 
the warp dDrim mel & Spergel , 2001 ; Lopez-Corredoira et al. 



2002t iMomanv et al.L I2006i) . and on the other hand the scale- 
height of the stars can show differences up to 50 per cent of 
the value for h^Rp,) in the range 5 < < 10 kpc dAlyes . 
2OOOI; iLopez-Corredoira et al.L l2002t iMomany et"all l2006h to 
fit a global Galactic disc model which accounts for all these in- 
homogeneities is, at least, tricky. It is because of this than the 
results in the Galactic model parameters might depend on the 
sample of Galactic coordinates used, as the combined effect of 
the warp and flare will be different at different directions in the 
Galaxy, hence at different lines of sight. 

In this paper, we derived the structural parameters of the 
thin and thick discs of the Milky Way by using data from 
2MASS survey, to observe changes in the parameters with the 
Galactic longitude. We briefly describe the 2MASS data, den- 
sity functions and the estimation of the Galactic model parame- 
ters in Section 2. Dependence of the Galactic model parameters 
with the Galactic longitude is given in Section 3. Finally, our 
main results are discussed and summarized in Sections 4 and 
5, respectively. 

2. 2MASS data 



The Two Micron All Sky Survey (2MA 55. 1 Skrutskie et al.Ll2006h 



provides the most complete database of near infrared (NIR) 
Galactic point sources available to date. During the develop- 
ment of this survey two highly-automated 1.3-m telescopes 
were used, one at Mt. Hopkins (AZ, USA) to observe the 



In ILopez-Corredoira et alJ (l2002h was presented a method of 
deriving stellar densities and the interstellar extinction along a 
given line of si ght, which was developed in a subsequent series 
of p apers J Drimmel. Cabrera-Lavers & Lopez-Corredoir; 
"20031 



Picaud. Cabrera -Lavers & Garzon . 



doirai 
20031: 



Lopez-Corredoira et al.. 2004). The method relies in extracting 
the red clump (RC) population from the colour-magnitude 
diagrams (CMDs) as they are the domi nant giant population 



dCohen et ak. .2000. ; Hammerslev et al. , I2OOO.) . These stars 



can be easily isolated as they form a conspicuous feature 
in the CMDs. The absolute magnitude and intrinsic colour 
of this population are weU defined: M(K) - -1.62 + 0.03, 
(J - K)o = 0.61 + 0.01 w it h a s r nall dependence with the 
metall icity and age (Alves, '2000'; 'Grocholski & Sarajedinij 



2002; ISalaris & Girardi . 2002; Pietr zyrisk i et al .. 2003 ). Of 



course, there is some dispersion with respect to those values 
(cr[(7 - K)()] ~ 0.1), but we know that this population is 
dominant and the dispersion of absolute magnitudes and 
colours is not very large so the use of an average value of 
M{K) and (J - K)() is a good approximation. Therefore, we 
can extract spatial information directly from the apparent 
magnitudes and (J - K) colours of the RC stars in the CMDs. 



The method has been used before in Cabrera-Lavers et al 



d2005h to analyze the thick disc of the Milky Way, extrac- 
ting some results about the structural parameters of this com- 
ponent that were in good agreement with previous estimates 
for the Galactic thick disc. The full details of the applica- 
tion of the method can be fo und in that paper and also in 
Lopez-Corredoira et alj d2002l) so they will not be repeated 



here. As a brief summary, the RC stars are firstly isolated in 
the CMDs by means of theoretical tr aces predicted using the 
"SKY" model dWainscoat et al. . 1992 ). that define the approx- 
imate area in the CMD where the RC population lies. Once 
these stars are identified, RC stars are extracted around a trace 
fitted to the maxima of a series of consecutive Gaussian fits to 
different colour histograms at fixed apparent magnitudes, niic 
(see Fig. [TJ. The distance along the line of sight to those stars 
can be derived easily from the apparent colour and magnitude. 
By obtaining the number of RC stars in each interval of appar- 
ent magnitude, the stellar density can be derived once the abso- 
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lute magnitude and intrinsic colour of the RC are known (that 
is the main assumption of the method). As we obtain densities 
along the line of sight, they can be transformed into densities 
in cylindrical coordinates (R, z). 

-5°<l<5° 60°<b<70° 



10 
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Fig. 1. K vs. J - K 2MASS colour-magnitude diagram for the 
field centred at Z = 0°. The fitted trace that we assign to 
RC giants population (solid line) and the limits for extrac- 
tion of the stars (dashed lin es) are also shown (see Sect. 2 in 
Cabrera-Lavers et alj (120051) for details). 



The possible contamination due to dwarf stars in the ex- 
tracted counts is the more restrictive issue when applying the 
method. To minimize this effect, we extracted only stars up 
to niK <10 from the CMDs. In this range, contamination is 
less than 10 per cent of the total numb er stars, and it is even 



lower for brighter apparent magnitudes (iCabrera-Lavers et al. 



20051) . A magnitude mk = 10 corresponds to a distance from 
the Sun of approximately 2.5 kpc for a RC star, suitable for a- 
nalyzing the disc more than 1 kpc above the Galactic plane but 
unable to reach the range of distances where the halo dominates 
the counts. Thus, the following analysis is made ignoring this 
Galactic component. 

Other aspects that affect the method, as the metallicity de- 
pendence of the colour of the RC population, any possible 
contamination o f lower luminosi ty giants in the counts, or the 
Malmquist bias jMalmquistl [l92 0) in the absolute magnitude 
of the RC stars have been taken into account, and their pos- 
sible effects in the re s ults w ere deeply discussed in Sect. 3 of 
Cabrera-Lavers et alJ (l2005h . In any case, the uncertainties are 
always lower than the one coming from the contamination of 
dwarf stars. 

We have first collected 10x10 degree fields in the 2MASS 
catalogue centered at fixed Galactic longitudes (in Al = 10° 
bins), with 60° < b < 70°. We have to average several fields to 
obtain a sufficient number of RC stars to build up a conspicuous 
feature in the CMD which could be easily identifiable. The den- 
sities extracted are fitted by a double exponential following eq. 



([T]i, for both the thin and thick discs, with the additional con 
straint of producing the same local disc space density for the 



thin disc as it was obtained previously in lCabrera-Lavers et al 
(2005). 

Therefore, we assume for the discs a distribution as follows: 

Pi(R, z) = n, exp(-\z\/h,j) exp(-(R - Ro)/Hi), (1) 

where z - Ze + rsm{b), r is the distance to the object 
from the Sun, b is the Galactic latitude, Zq is the vertical dis- 
tance of the Sun from the Galactic plane we assume to be 
15 pc jHammerslev et al. , 1 1994 ). R is the projected of the 
Galactocentric distance on the Galactic plane, Rq is the solar 
distance from the Galactic center (Reid, 1993, 8 kpc), /z, , and 
Hi are the scaleheight and scalelength, respectively, and n, is 
the normalized density at the solar radius. The suffix / takes the 
values 1 and 2 as long as the thin and thick discs are considered. 

As this study focuses on the dependence of the scale- 
height and solar normalization on the Galactic longitude, for 
the scalelengths of the thin and thick discs we used the val- 
ues of 2.1 kpc dLopez-Corredoira et al 1 i2()02l) and 3 kpc 



(ICabrera-Lavers et ; 



■^opez-Corredoira et al.L l/Ul)2l) and i kpc 
al.l 2005h respectivelv. which were also ob 



tained with 2MASS data. In any case, it has to be noted that 
the range of Galactocentric distances covered along each line 
of sight is so small that the effec t of changing the value of th e 



radial scalelengths is negligible (ICabrera-Lavers et al.L l2005h . 
The possible implications of a different scalelength of the thick 
disc will be discussed in Section l4r2l 

Fig. |2] shows stellar densities profiles extracted at four dif- 
ferent Galactic longitudes, as an example of the application of 
the method, while the Galactic model parameters are given in 
Table [T] Columns (2) and (3) give the scaleheight of the thin 
and thick discs, respectively, while the local space density of 
the thick disc («2/ni) is given in Col. (4). Error bars in the 
scaleheights come from the uncertainty in the double exponen- 
tial fit, by assuming Poisson noise in the extracted counts and 
taking into account the uncertainty in the distance above the 
plane as we have grouped 10x10 degree fields. 

3. Dependence of the Galactic model parameters 
with the Galactic longitude 

Fig. [3] shows the variation of the scaleheight of the thin 
disc with the Galactic longitude. No obvious global trend 
is observed in the thin disc scaleheights, with a sam- 
ple of values well compatible with a constant value of < 
h, >=187 (cr=36) pc, well in agreement with us u al es- 



timates for thi s component (e.g.. [ Bah call & Soneir al 11984 



iRob in & Crezl Il986l; iReid & Maiewski .1993.: Siegel et al 
2002; Lopez-Corredoira et all 2002h . As shown in Fig. [3] 



have overplotted a sinusoidal fit, no physical meaning, to the 
data. One can say that the scaleheight of the thin disc shows 
slight variation with the Galactic longitude. However, the scat- 
ter and the error bars are large, being compatible with a con- 
stant value independent of the Galactic longitude. Hence, we 
can conclude that the thin disc scaleheight is found to be a con- 
stant as a function of Galactic longitude. 

For the thick disc the variation of the scaleheight with 
the Galactic longitude is more pronounced (Fig. |4li. There is 



4 



Cabrera-Lavers et al.: Galactic model parameters 



-5"<l<5" 60"<b<70" 



-8.5 
-9 
-9.5 
-10 
-10.5 
-11 
-11.5 
-12 
-12.5 
-13 
-13.5j 

-8.5 
-9 
-9.5 
-10 
-10.5 
-11 
-11.5 
-12 
-12.5 
-13 
-13.5, 





-8.5 r 




-9 




-9.5 




-10 




-10.5 


a. 




c 


-11 


_i 






-11.5 




-12 




-12.5 




-13 




-13.5^ 







-8.5 
-9 
-9.5 
-10 
-10.5 
-11 
-11.5 
-12 
-12.5 
-13 
-13.5, 



Table 1. Galactic model parameters for the 36 star fields ana- 
lyzed with 2MASS data. 
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Fig. 2. Density of red clump stars (in stars pc ^) obtained from 
the CMDs for the Galactic longitudes of 0°, 90°, 170°, 270°. 
The best fits for a sum of thick disc (dashed line) and thin disc 
(dotted line) components are also shown as a solid hne. Error 
bars come from the Poisson noise in the extracted counts and 



Fig. 3. Variation of the scaleheight of the thin disc with the 
Galactic longitude. 
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a maximum in the inner Galaxy (|/| < 30°), whereas in be- 
tween 150° < I < 230° there is a local minimum. The 
mean of the scaleheight for the data, < h~ >- 957 (cr=100) 
pc is well w i thin t he given range fo und in the literature 
(ISpagna et all 1 19961: fNg et all 1 19971: iBuser. Rons & Karaali . 



< b > = 50 



1998'. '1999^. 'Siegel et al.'. ' 20021: iLarsen & Humphreysl 12003 
Cabrera-Lavers et aL.2005i) . 
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Fig. 4. Variation of the scaleheight of the thick disc with the 
Galactic longitude. 



The variation of the local space density of the thick disc 
relative to the local space density of the thin disc (na/ni) with 
the Galactic longitude is given in Fig.|5] The trend gives a con- 
stant local space density, (n2/«i) ~ 10 per cent, which confirms 
that the variation in h- in Fig.|4]is not an artifact of an anticor- 
relation in the local normalization (assuming, for example, a 
constant value (ni/ni) = 10% we still get the result in Fig.|4]i. 

We checked the effect of the Galactic latitude on the varia- 
tion of the scaleheight of the thick disc with the Galactic longi- 
tude for the 2MASS data. Fig.|6]shows the variation of the scale- 
height with the Galactic longitude for fields at< b >- 50°. The 
trends in Figs. |4] and |6] are similar, although the zero points are 
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Fig. 5. Variation of the relative space density of the thick disc 
with the Galactic longitude. 
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Fig. 6. Variation of the scaleheight of the thick disc with the 
Galactic longitude, for fields with Galactic latitude < b >- 50°. 



different. The equations of the sinusoidal fits for the data in the 
fields with < b >- 50° and < b >- 65° are as follows: 

h- = 918 - 141.45 sin (/ - 84.65°), (at < b 50°), (2) 

/!, = 959- 108.02 sin (Z- 63.87°), (at < b >^ 65°). (3) 

The difference between the scaleheights is larger (it amounts 
up to ~ 100 pc) for the longitudes corresponding to different 
minima in the two figures. No evident differences are found 
between the overall shape of the distributions. Hence the ob- 
served behaviour is associated to a global trend in the thick 
disc rather than to something exclusive of a restricted range of 
latitudes. 

4. Discussion 

4.1. Does the flare affect the scaleheight? 

It is well known that the Galactic disc shows a flare, which 
produces an incr ease in the scaleheight a s we rn ove outwards 
in the Galaxy. In iLopez-Corredoira et al.l (120021) the flare was 
modelled by an exponential increase of with R, obtaining a 
law that approximately reproduced the observed counts up to 
R < 15 kpc. On the contrary, for the inner Galaxy a flare with 
the opposite trend (an increase in h, as R decreases) was found 
by lLopez-Corredoira et al. I (l2004 . who then proposed an ex- 
pression that summarizes both regimes with a smooth transition 
between them: 



h,(R) = h,(Re) X [1 -H 0.21(7? - /Jo) + 0.056(7? - RqY] 



(4) 



The variation of the scaleheight of the thin and thick disc 
with Galactic longitude could be related with the flare itself. 
When the line of sight is pointing to the inner Galaxy (\l\ < 30°) 
the mean Galactocentric distance of the sources is lower than 
when the line of sight is in the anticentre direction. We can 
then translate the h, vs. longitude plot in a /i- vs. R plot by 
deriving the mean Galactocentric distance to either the thin or 
thick disc sources, a procedure that is made by binning the data 
by the Galactocentric distance, R, and determining the mean 
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Fig. 7. Variation in the scaleheight of the thin (above) and thick 
(below) discs with R. 



radii of the red clump stars in each bin to average them. Figure 
|7] shows the variation in the scaleheight of the thin and thick 
disc with Galactocentric distance. It has to be noted that the 
range of Galactocentric distances is very small, so a possible 
derivation of parameters for the flare is far from being conclu- 
sive. We have compared Lopez-Corredoira et al.'s law in both 
cases (dashed lines in the plot) considering /i; (Rq) = 200 pc 
and h~ (Rq) - 900 pc for the thin and thick disc, respectively. 
While the predictions for the thin disc are compatible with the 
data (although a constant value of h- can also fit the data), for 
the thick disc the observed trend is just the opposite with an 
increase of when moving to lower values of R (represented 
by means of a linear fit on the lower panel of Fig. IT), a result 
that is in complete disagreement with the obtained for the outer 
thin disc (R>6 kpc) in Lopez-Corredoira et al. (2002). 

To obtain a more accurate representation of the flare in 
the thick disc it would be necessary to increase the range 
of Galactocentric distances. For this reason, 2MASS data at 
< b >= 50° are useful. They correspond to lines of sight at 
lower heights above the plane, thus the mean Galactocentric 
distance of the thick disc sources increases improving the fit. 
In Fig. [8] we show the variation in the scaleheight of the thick 
disc with respect to the mean Galactocentric distance, com- 
bining the results at < b >- 50° (filled circles) with those at 
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Fig. 8. Variation of the scaleheight of the thick disc with R by 
using 2MASS data at < 50° and < >= 65°. The best 

linear fit to the observed trend is also shown. 



< b >- 65° (open circles). By means of a linear fit to the data 
we obtain the following expression: 

h,(R) ^(940 + 20) X [I- (0.12 ±Om)(R-RQ)] (pc) (5) 

It is not so common to observe an increase of the scale- 
height as we move to the inner Galaxy. Furthermore, no simi- 
lar analysis of a possible flare in t he thick disc is found in the 
literature. In a very recent article. iMomanv et al.l (120061) anal- 
ysed the warp and flare of the disc of the Milky Way but only 
concentrated in those related to the thin disc. By using red gi- 
ant branch stars in the range 0.45 < z < 2.25 kpc they ob- 
tained a slight increase of the scaleheight of the disc with in- 
creasing Galactocentric radii, up to R ~15 kpc (their fig. 15). 
They claimed that the overall mean scaleheight around 0.65 
kpc they obtained was a reflection of a mixture of the thin and 
thick disc populations. If we put together the known trend in 
the scaleheight of the thin disc (that is, an increasing /i- as R 
increases) with the result obtained in this paper for the thick 
disc (an increase in h~ as R decreases), one must expect that a 
sort of 'cancellation' appears, producing a nearly constant h, at 
intermediate heights above the plane. 

To check this possibility, we have made a simple simu- 
lat ion for a combination of a thi n disc with the density laws 
of iLopez-Corredoira et al. I (l2004 and a thick disc with an in- 
crease in represented by eq. (|5]l and a relative space density 



of 10 p er cent of the thin disc according to lCabrera-Lavers et al 
(12005b . The resulting densities are fitted by a single expo- 



nential law (as if only one component were present), obtain- 
ing in that way the expected variation in scaleheight with the 
Galactocentric distance. 

When com paring this simulation with the data taken from 
Momanv et al.l (2006) (Fig. |9]l a good agreement is observed. 
In the range of distances 4 < R <9 kpc a nearly constant scale- 
height is expected that is not far from that observed (note we 
did not perform any kind of fit to Momany et al.'s data, so the 
coincidence in the value for h~at4<R<7 kpc is very interes- 
ting). Predicted values for R > 10 kpc are larger than those, but 
it was expected due to the excessively high increase in for 
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the thin disc that Lopez-Corredoira et al.'s model produces in 
this range of Galactocentric distances. A more modest increase 
in h, with R for the th in disc (as obtained in Alves (2000) or in 
Momany et al. liooi)) would produce a better agreement with 



the data, but this was beyond the scope of this comparison. For 
R < 4 kpc, again an excessive increase in is obtained in the 
simulation, but again it was known that the validity of Lopez- 
Corredoira et al.'s model for R < 2.5 kpc is not well demon- 
strated. 

We have also introduced some changes in the simulation. 
By eliminating the contribution of the thick disc (that is, as- 
suming that only the thin disc is observed) we recovered the 
predictions of the Lopez-Corredoira et al.'s model, with the 
scaleheight corresponding to the thin disc (dashed line in Fig. 
|9]l. It is even more remarkable how when a flaring in the thick 
disc is neglected (assuming a constant scaleheight of 940 pc for 
this component) the trend obtained still follows the increase in 
the scaleheight due the thin disc, but with higher h, values as 
corresponds with the mixing of both populations (grey solid 
line in Fig.|9]l. The result is also less sensitive to other parame- 
ters, as the relative space density of the thick disc (no changes 
were obtained by varying it from 5 to 15 per cent) or the value 
of Rq (assumed as 8 kpc). 

A combination of a flared thin disc and a flared thick disc 
with opposite trends seems, then, to be compatible with the 
observed scaleheights at intermediate heights above the plane, 
as those of Momany et al. (2006). Thus, the observed variation 
in the scaleheight of the thick disc with the galactic longitude 
has been obtained here as expected. 

4.2. Effect of the scalelength of the thick disc 

In Section 12.11 we addressed what we considered a fixed 
scalelength of 3 kpc for the thick disc, as obtained in 



Cabrera-Lavers et al.l ( 120051) . However, a different scalelength 
would reproduce the observed variation in the scaleheight with- 
out the need for a flare in the disc itself. For example, a shorter 
scalelength of 2 kpc for the thick disc would increase the pre- 
dicted star counts in the inner disc and reduce them in the outer 
disc, relative to the assumed 3 kpc scalelength, simulating a 
higher scaleheight in the inner disc and a lower scaleheight in 
the outer disc, as was found in Section 4. 

However, the effect of the scalelength on the derived den- 
sities is negligible, as we are moving in a very narrow range 
of Galactocentric distances. To reproduce the same density as 
that resulting from a scalelength of 3 kpc and a scaleheight of 
11 00 pc, like those derived in the inner Galaxy but with a con- 
stant scaleheight of 800 pc, the scalelength has to change by 
a factor ~1.4 with respect to the assumed value of 3 kpc (that 
is, around 4.2 kpc), a value that is 11 cr above that derived in 
ICabrera-Lavers et all (l2005b . 

In order to check this statement, we repeated the analysis 
of Section 2 in a series of nine test fields but now assuming 
two different scalelengths for the thick disc: H - 2 kpc and 
H = 2.5 kpc, respectively. The results are shown in Fig. [TOl 
together with sinusoidal fits to the data. The observed trend 
is nearly coincident with that of Fig. |4]in both amplitude and 
the positions of the minima. Hence, the effect of a different 
scalelength from the assumed one of 3 kpc has to be discarded 
as being responsible for the results discussed in this paper. 



5. Conclusion 

By using NIR data from 2MASS we have observed significant 
changes in the structural parameters of the thin and thick discs 
of the Milky Way, that are more than simple fluctuations caused 
by model fitting to the data. 

A plausible explanation for the observed changes in the 
scaleheight of the discs arises from the combined effect of 
the warp and flare in the Galactic thin and thick disc, which 
displaces both structures from the assumed position expected 
for a "smooth disc" (that is, a flat, azimuthally homogeneus 
and constant scaleheight disc). This implies that in some di- 
rections the effect is more noticeable than in others, produc- 
ing variations in the parameters resulting from a fit of a sin- 
gle Galactic model. The observed trend in the change of the 
scaleheight with galactic longitude is well reproducible the- 
oretically by assuming the accretion of intergalactic matter 
on to the Galactic disc, which in general produces a differ- 
ent pressure depending on both the Galactocentric radius and 
azimuth, the same pressure which would also be responsi- 
ble for the formation of S-warps and/or U-war p s in g alaxies 
( Lopez-Corredoira. Betancort-Rijo & BeckmaH 20021) . The 
pressure would be similar to a piston mechanism, only from 
one side of the disc ( Sanchez-Salcedoi I2OO6, §4.6). As a result 
of this mechanism the scaleheight of the disc would present 
a similar shape to that obtained in Figure |4] However, it is far 
from our intention to consider that our result is a demonstration 
that this speculative hypothesis is the correct one in explaining 
the uncertain origin for the Milky Way's flare and warp; but it 
is consistent with the observed data. 
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Fig. 10. Variation in the scaleheight of the thick disc with the 
Galactic longitude obtained by assuming a scalelength of H-2 
kpc (above) or H-2.5 kpc for this component (below). 



The flare obtained in the thick disc from 2MASS data has 
an opposite trend to that commonly assumed for the thin disc, 
with a scaleheight that increases as we move to the innermost 
Galaxy, being the first time that a possible flare in this Galactic 
component is observed. It is beyond the scope of this paper 
to determine from the observed results a theoretical scenario 
that supports the origin of this possible flare, although some 
asymmetr ies in the thick dis c hve been well known for some 



time now (IParker et aUl2003h . We suspect that the different be- 



haviour in the inner Ga laxy is related to the interacti on of the 4 



kpc inner Galactic bar dLopez-Corredoira et all 12006) with its 



surroundings. However, it is more difficult to explain the pos- 
sible effect of this component on the thick disc, although some 
known asymmetries of this component ha ve been sugg e sted a s 



having been caused by the Galactic bar dParker et al.L |2004|) 



Unfortunately, we do not have enough information to add any- 
thing in this regard. Probably, kinematic analyses by using ra- 
dial velocity data as those provided by the Sloan Extension for 
Galactic Understanding and Exploration (SEGUE, Newberg et 
al., 2003) or the Radial Velocity Experiment (RAVE, Steinmetz 
et al., 2006) will be very valuable to go deeper in this topic. 
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